Stable isotope composition of pesticides in commercial formulations: The ISOTOPEST database.

By assessing the changes in stable isotope compositions within individual pesticide molecules, Compound Specific Isotope Analysis (CSIA) holds the potential to identify and differentiate sources and quantify pesticide degradation in the environment. However, the environmental application of pesticide CSIA is limited by the general lack of knowledge regarding the initial isotopic composition of active substances in commercially available formulations used by farmers. To address this limitation, we established a database aimed at cataloguing and disseminating isotopic signatures in commercial formulations to expand the use of pesticide CSIA. Our study involved the collection of 25 analytical standards and 120 commercial pesticide formulations from 23 manufacturers. Subsequently, 59 commercial formulations and 25 standards were extracted, and each of their active substance was analyzed for both δ13C (n = 84) and δ15N CSIA (n = 43). The extraction of pesticides did not cause significant isotope fractionation (Δ13C and Δ15N < 1‰). Incorporating existing literature data, stable carbon and nitrogen isotope signatures varied in a relatively narrow range among pesticide formulations for different pesticides (Δ13C and Δ15N < 10‰) and within different formulations for a single substance (Δ13C and Δ15N < 2‰). Overall, this suggests that pesticide CSIA is more suited for identifying pesticide transformation processes rather than differentiating pesticide sources. Moreover, an inter-laboratory comparison showed similar δ13C (Δ13C ≤ 1.2 ‰) for the targeted substances albeit varying GC-IRMS instruments. Insignificant carbon isotopic fractionation (Δ13C < 0.5‰) was observed after 4 years of storing the same pesticide formulations, confirming their viability for long-term storage at 4 °C and future inter-laboratory comparison exercises. Altogether, the ISOTOPEST database, in open access for public use and additional contributions, marks a significant advancement in establishing an environmentally relevant pesticide CSIA approach.

A Novel Multi-ion Evaluation Scheme to Determine Stable Chlorine Isotope Ratios (37Cl/35Cl) of Chlordecone by LC-QTOF.

Organochlorinated pesticides are highly persistent organic pollutants having important adverse effects in the environment. To study their fate, compound-specific isotope analysis (CSIA) may be used to investigate their degradation pathways and mechanisms but is currently limited to 13C isotope ratios. The assessment of 37Cl isotope ratios from mass spectra is complicated by the large number of isotopologues of polychlorinated compounds. For method development, chlordecone (C10Cl10O2H2; hydrate form), an organochlorine insecticide that led to severe contamination of soils and aquatic ecosystems of the French West Indies, was taken as a model analyte. Chlorine isotope analysis of chlordecone hydrate was evaluated using high-resolution liquid chromatography quadrupole time-of-flight mass spectrometry (LC-QTOF-MS), enabling smooth ionization to detect the molecular ion. First, a new evaluation scheme is presented to correct for multiple isotope presence in polychlorinated compounds. The scheme is based on probability calculations of the most frequent isotopologues, distributions by binomial probability functions, and corrections for the presence of nonchlorine heavy isotopes. Second, mobile-phase modifiers, ionization energy (sampling cone tension) and scan time were optimized for accurate chlorine isotope ratios. Chlordecone standard samples were measured up to 10-fold and bracketed with a second chlordecone external standard. δ37Cl values were obtained after conversion to the SMOC scale by a two-point calibration. The robustness of the analysis method and evaluation scheme were tested and gave satisfactory results with standard errors (σm) of ±0.34‰ for precision and ±0.89‰ for long-term accuracy of chlorine isotope ratios of chlordecone hydrate. This work opens perspectives for applications of the C-Cl CSIA approach to investigate the fate of highly toxic and low reactive polychlorinated compounds in the environment.

Chlorinated solvents source identification by nonlinear optimization method.

In this work, chloride ions were used as conservative tracers and supplemented with conservative amounts of chloroethenes (PCE, TCE, Cis-DCE, 1,1-DCE), chloroethanes (1,1,1-TCA, 1,1-DCA), and the carbon isotope ratios of certain compounds, the most representative on the sites studied, which is a novelty compared to the optimization methods developed in the scientific literature so far. A location of the potential missing sources is then proposed in view of the balances of the calculated mixing fractions. A test of the influence of measurement errors on the results shows that the uncertainties in the calculation of the mixture fractions are less than 11%, indicating that the source identification method developed is a robust tool for identifying sources of chlorinated solvents in groundwater.

Assessment of the contamination by 2,4,6-tribromophenol of marine waters and organisms exposed to chlorination discharges.

2,4,6-tribromophenol (TBP) is implied in the production of brominated flame retardants but is also a major chlorination by-product in seawater. A growing number of studies indicate that TBP is highly toxic to the marine biota, but the contribution of anthropogenic sources among natural production is still under question concerning its bioaccumulation in marine organisms. Here, several water sampling campaigns were carried out in the industrialized Gulf of Fos (northwestern Mediterranean Sea, France) and clearly showed the predominant incidence of industrial chlorination discharges on the TBP levels in water, at the 1-10 ng L-1 level in average and reaching up to 580 ng L-1 near the outlets. The bioaccumulation of TBP was measured in 90 biota samples from the Gulf of Fos. The concentrations found in European conger muscle tissues (140-1000 ng g-1 lipid weight, in average), purple sea urchin gonads (830-880 ng g-1 lipid weight, in average), and Mediterranean mussel body (1500-2000 ng g-1 lipid weight, in average) were above all published references. Significant correlations with fish length (European conger) and gonad somatic index (purple sea urchin) were also identified. Comparatively, fish, urchins and mussels from other Mediterranean sites analyzed within this study showed a lower bioaccumulation level of TBP, consistently with what found elsewhere. Industrial outflows were thus identified as hotspots for TBP in seawater and marine organisms. The environmental risk assessment indicated a high potential toxicity in the industrial Gulf of Fos, in particular near the outlets, and a limited threat to human but toxicological references are lacking.

Underestimation of Anthropogenic Bromoform Released into the Environment?

Bromoform (CHBr3) belongs to very-short-lived substances (VSLSs), which are important precursors of reactive bromine species (BrOx) contributing to tropospheric and stratospheric chemistry. To date, most models calculating bromine product emissions to the atmosphere only consider the natural production of CHBr3 from marine organisms such as macroalgae and phytoplankton. However, CHBr3 has many other anthropogenic sources (coastal industrial sites, desalination and wastewater plants, ballast waters, and seawater toilets) that may drastically increase the amounts emitted in the atmosphere. Here, we report the levels of CHBr3 released in water and air (according to real-time and offline measurements by proton-transfer-reaction time-of-flight mass spectrometry (PTR-ToF-MS) and gas chromatography with electron capture detection (GC-ECD)) in a highly industrialized area where 3 million cubic meters of chlorinated seawater is released each day, which were measured during six field campaigns (at sea and on land) distributed over 3 years. The highest levels found during this survey (which were correlated to the physical-chemical characteristics of the water, meteorological and hydrological conditions, salinity, and temperature gradients along the water column) reached 34.6 µg L-1 in water (100-10 000 times higher than reported natural levels) and 3.9 ppbv in the air (100 times higher than the maximum reported value to date). These findings suggest the need to undertake sampling and analysis campaigns as close as possible to chlorinated discharges, as anthropogenic CHBr3 sources from industrial discharges may be a missing factor in global flux estimates or organic bromine to the atmosphere.

Using 1,1,1-Trichloroethane degradation data to understand NAPL dissolution and solute transport at real sites.

Analytical and numerical models describing the evolution of contaminant concentrations in the plume associated with the dissolution of NAPL source and degradation processes were presented in the literature. At real sites and particularly in complex aquifers like chalk, it is difficult to understand how the sources of contaminants evolve with time. 1,1,1-Trichloroethane (1,1,1-TCA) is one of the few compounds with a well-known hydrolysis constant, that can help to improve knowledge of the contaminant sources and transport rates of dissolved contaminants in groundwater by dating the spill. In this work, different scenarios that could explain the evolution of the concentrations of 1,1,1-TCA and its degradation product 1,1-Dichloroethene (1,1-DCE) at a real contaminated site were investigated by analytical and numerical modelling. The results show that (1) the peaks of concentration time series do not correspond to a single contamination event even in the case of a complex medium, (2) the multiphasic behavior of the concentration time series is dictated by the dissolution in a heterogeneous medium, and (3) the persistence of the concentrations can arise from a small residual organic phase or transport in dual domain medium.

Passive sampling of chlorinated paraffins by silicone: Focus on diffusion and silicone-water partition coefficients.

Short-chain chlorinated paraffins (SCCPs) are under regulation through the European Water Framework Directive and were recently classified as POPs. Consecutively, the increasing use of middle-chain chlorinated paraffins (MCCPs) becomes of growing concern. Knowledge on the occurrence of chlorinated paraffins (CPs) is still scarce particularly in water phase. To achieve sufficient method sensitivity, the passive sampling approach, acting as a relevant alternative to usual grab sampling, has been considered only very recently for the monitoring of CPs in water. The present work aimed at determining the diffusion coefficients in silicone (Ds) and the silicone-water partition coefficients (Ksw) of various CP groups, having different chlorine contents and carbon chain lengths, in four commercial CP mixtures. Log Ds (-10.78 to -10.21) was found to vary little and to be high for the groups of CPs studied. Thus, their uptake in silicone is controlled by the water boundary layer, which allows to consider the release of performance and reference compounds for in-field estimation of the sampling rate. Moreover, CPs partitioned strongly towards silicone rubbers. Both the chlorination degree and the carbon chain length of CPs cause large uncertainties in the partitioning between silicone and water (log Ksw between 4.85 and 6.30), indicating that instead of an average value, differentiated Ksw should be used to estimate aqueous CPs more accurately. Even so, the probable influence of chlorine atoms position on polarity and partitioning may be an argument for favoring sampling in the kinetic stage.

Dehalogenation of α-hexachlorocyclohexane by iron sulfide nanoparticles: Study of reaction mechanism with stable carbon isotopes and pH variations.

The biodegradation of hexachlorocyclohexanes (HCHs) is known to be accompanied by isotope fractionation of carbon (13C/12C), but no systematic studies were performed on abiotic degradation of HCH isomers by iron (II) minerals. In this study, we explored the carbon isotope fractionation of α-HCH during dechlorination by FeS nanoparticles at different pH values. The results of three different experiments showed that the apparent rate constants during dehalogenation of α-HCH by FeS increased with pH. The lowest apparent rate constant value α-HCH during dehalogenation by FeS was 0.009 d-1 at pH value of 2.4, while the highest was 1.098 d-1 at pH 11.8. α-HCH was completely dechlorinated by FeS only at pH values 9.9 and 11.8, while the corresponding apparent rate constants were 0.253 d-1 and 1.098 d-1, respectively. Regardless of the pH used, the 1,2,4-trichlorobenzene (1,2,4-TCB), 1,2-dichlorobenzene (1,2-DCB), and benzene were the dominant degradation products of α-HCH. An enrichment factor (εC) of -4.7 ± 1.3‰ was obtained for α-HCH using Rayleigh model, which is equivalent to an apparent kinetic isotope effect (AKIEC) value of 1.029 ± 0.008 for dehydrohalogenation, and of 1.014 ± 0.004 for dihaloelimination, respectively. The magnitude of isotope fractionation from this study suggests that abiotic isotope fractionation by FeS must be taken into account in anoxic sediments and aquifers contaminated with HCH isomers, when high concentrations of FeS are present in the above-mentioned anoxic environments.

Effect of NAPL mixture and alteration on 222Rn partitioning coefficients: Implications for NAPL subsurface contamination quantification.

Soils and groundwater are often contaminated by complex organic mixtures also called Non Aqueous Phase Liquids (NAPLs). Several techniques such as drilling, monitoring of soil gas or injection of tracers are traditionally used to quantify NAPLs in aquifers but are complex to perform. The use of natural soil gas such as 222Rn could be an easy and cheap alternative. This method requires the knowledge of the radon NAPL-water partitioning coefficients (Kn-w.). Once spilled on soil, NAPL will undergo degradation (evaporation, effects of sun light among others) and this degradation could impact the Kn-w. This study aims at investigating the partitioning coefficients of complex NAPLs such as commercial diesel fuel and gasoline in relation to degradation such as evaporation and UV-degradation. For that purpose, batch experiments and GCMS investigations were carried out. The results show different Kn-w for the commercial diesel fuel (60.7 ± 6.1) and gasoline (37.4 ± 5.6). The results also show different Kn-w behaviors in relation with degradation. Degraded diesel fuel display opposite Kn-w values (74.8 ± 7.5 and 25.1 ± 2.5 for UV degraded and evaporated diesel fuel, respectively), compared to fresh one. Degraded gasoline shows no significant variations of the Kn-w compared to fresh one. The molecular investigation reveals the removal of the most volatile fraction for the evaporation treatment, whereas UV-degradation do not have pronounced effects on the chromatogram pattern. For the gasoline molecular investigation, no difference is observed between the treatments excepted a very slight removal of the lightest compounds under evaporation. These results show that NAPL degradation have effects on the Kn-w for diesel fuel and no significant effects for gasoline, at least with these degradation paths. This Kn-w variation will have in fine effects on 222Rn activity interpretation and NAPL subsurface quantification.

Change in Evaluation Mode Can Cause a Cheerleader Effect.

The cheerleader effect describes the phenomenon whereby faces are perceived as being more attractive when flanked by other faces than when they are perceived in isolation. At least four theories predict the cheerleader effect. Two visual memory processes could cause a cheerleader effect. First, visual information will sometimes be averaged in the visual memory: the averaging of faces could increase the perceived attractiveness of all the faces flanked by other faces. Second, information will often be combined into a higher-order concept. This hierarchical encoding suggests that information processing causes faces to appear more attractive when flanked by highly attractive faces. Two further explanations posit that comparison processes cause the cheerleader effect. While contrast effects predict that a difference between the target face and the flanking faces causes the cheerleader effect due to comparison processes, a change in the evaluation mode, which alters the standard of comparison between joint and separate evaluation of faces, could be sufficient for producing a cheerleader effect. This leads to the prediction that even when there is no contrast between the attractiveness of the target face and the flanking faces, a cheerleader effect could occur. The results of one experiment support this prediction. The findings of this study have practical implications, such as for individuals who post selfies on social media. An individual's face will appear more attractive in a selfie taken with people of low attractiveness than in a selfie without other people, even when all the faces have equally low levels of attractiveness.

Intramolecular non-covalent isotope effects at natural abundance associated with the migration of paracetamol in solid matrices during liquid chromatography.

Position-specific isotope analysis by Nuclear Magnetic Resonance spectrometry was employed to study the 13C intramolecular isotopic fractionation associated with the migration of organic substrates through different stationary phases chromatography columns. Liquid chromatography is often used to isolate compounds prior to their isotope analysis and this purification step potentially alters the isotopic composition of target compounds introducing a bias in the later measured data. Moreover, results from liquid chromatography can yield the sorption parameters needed in reactive transport models that predict the transport and fate of organic contaminants to in the environment. The aim of this study was to use intramolecular isotope analysis to study both 13C and 15N isotope effects associated with the elution of paracetamol (acetaminophen) through different stationary phases and to compare them to effects observed previously for vanillin. Results showed very different intramolecular isotope fractionation profiles depending on the chemical structure of the stationary phase. The data also demonstrate that both the amplitude and the distribution of measured isotope effects depend on the nature of the non-covalent interactions involved in the migration process. Results provided by theoretical calculation performed during this study also confirmed the direct link between observed intramolecular isotope fractionation and the nature of involved intermolecular interactions. It is concluded that the nature of the stationary phase through which the substrate passes has a major impact on the intramolecular isotopic composition of organic compounds isolated by chromatography methods..

Quantification of Lambda (Λ) in multi-elemental compound-specific isotope analysis.

In multi-elemental compound-specific isotope analysis the lambda (Λ) value expresses the isotope shift of one element versus the isotope shift of a second element. In dual-isotope plots, the slope of the regression lines typical reveals the footprint of the underlying isotope effects allowing to distinguish degradation pathways of an organic contaminant molecule in the environment. While different conventions and fitting procedures are used in the literature to determine Λ, it remains unclear how they affect the magnitude of Λ. Here we generate synthetic data for benzene δ2H and δ13C with two enrichment factors εH and εC using the Rayleigh equation to examine how different conventions and linear fitting procedures yield distinct Λ. Fitting an error-free data set in a graph plotting the δ2H versus δ13C overestimates Λ by 0.225%⋅εH/εC, meaning that if εH/εCis larger than 22, Λ is overestimated by more than 5%. The correct fitting of Λ requires a natural logarithmic transformation of δ2H versus δ13C data. Using this transformation, the ordinary linear regression (OLR), the reduced major-axis (RMA) and the York methods find the correct Λ, even for large εH/εC. Fitting a dataset with synthetic data with typical random errors let to the same conclusion and positioned the suitability of each regression method. We conclude that fitting of non-transformed δ values should be discontinued. The validity of most previous Λ values is not compromised, although previously obtained Λ values for large εH/εC could be corrected using our error estimation to improve comparison.

Multistep partitioning causes significant stable carbon and hydrogen isotope effects during volatilization of toluene and propan-2-ol from unsaturated sandy aquifer sediment.

This study aimed at investigating whether stable isotopes can be used to monitor the progress of volatile organic compounds (VOCs) volatilization from contaminated sediment during venting. Batches of a dry aquifer sediment were packed into stainless steel HPLC columns, humidified with distilled water and later contaminated by either liquid toluene or propan-2-ol. The VOCs were then volatilized by a stream of gas at room temperature, and the concentrations and stable isotope ratios of gaseous VOCs were recorded by isotope-ratio mass spectrometry. During early stages of volatilization of toluene, the isotope ratios Δδ13C shifted to more negative values by about -3 to -5‰ and the Δδ2H by more than -40‰, while the concentration remained at or near initial saturated vapor concentration. Depletion of the isotope ratios in the gas was explained by the vapor-liquid fractionation process, which is amplified by successive self-partitioning steps of gaseous VOC into remaining liquid VOC. For propan-2-ol the carbon isotope shift was negative like for toluene, whereas the H shift was positive. Hydrogen bonding in the liquid propan-2-ol phase causes a normal vapor-liquid H isotope effect which was described already in classical literature. The isotope shifts in the present experiments are larger than previously reported shifts due to phase-change processes and reach the magnitude of shifts usually observed in kinetic isotope fractionation.

Intramolecular isotope effects during permanganate oxidation and acid hydrolysis of methyl tert-butyl ether.

Stable isotopes have been widely used to monitor remediation of environmental contaminants over the last decades. This approach gives a good mechanistic description of natural or assisted degradation of organic pollutants, such as methyl tert-butyl ether (MTBE). Since abiotic degradation seems to be the most promising assisted attenuation method, the isotopic fractionation associated with oxidation and hydrolysis processes need to be further investigated in order to understand better these processes and make their monitoring more efficient. In this study, position-specific isotope effects (PSIEs) associated with permanganate oxidation and acid hydrolysis of MTBE were determined using isotope ratio monitoring by 13C Nuclear Magnetic Resonance Spectrometry (irm-13C NMR) combined with isotope ratio monitoring by Mass Spectrometry (irm-MS). The use of this Position-Specific Isotopic Analysis (PSIA) method makes it possible to observe a specific normal isotope effect (IE) associated with each of these two abiotic degradation mechanisms. The present work demonstrates that the 13C isotope pattern of the main degradation product, tert-butyl alcohol (TBA), depends on the chemical reaction by which it is produced. Furthermore, this study also demonstrates that PSIA at natural abundance can give new insights into reaction mechanisms and that this methodology is very promising for the future of modeling the remediation of organic contaminants.

Laboratory-scale experimental and modelling investigations of 222Rn profiles in chemically heterogeneous LNAPL contaminated vadose zones.

The potential of LNAPL delineation by 222Rn soil-gas monitoring in a chemically heterogeneous vadose zone was investigated in this study based on laboratory (batch and columns) experiments and numerical modelling. An enhanced version of the MIN3P reactive transport code was used to simulate Rn transport in both uncontaminated and NAPL-contaminated vadose zones and results were validated against analytical solutions and laboratory experiments. Results show that 222Rn activity profiles are mainly controlled by porous media 222Rn production, vadose zone fluid saturations and especially the type and distribution of NAPL in contaminated areas. The results also show that decreases in 222Rn activity and variations in activity gradients provide evidence for the presence and saturation of NAPL. This study demonstrates that LNAPL delineation via 222Rn gas surveys at contaminated sites works best, if gas measurements extend as deep as possible and include regions where 222Rn activity decreases due to elevated NAPL content. In addition, collection and analysis of depth-discrete gas samples allows the characterization of vertical NAPL distribution based on the 222Rn activity gradient. The determination of 222Rn production in the unsaturated zone, as well as water capillary pressure curves are of key importance in enabling the discrimination of an uncontaminated from a NAPL-contaminated area.

Expanded uncertainty associated with determination of isotope enrichment factors: Comparison of two point calculation and Rayleigh-plot.

The enrichment factor (ε) is a common way to express Isotope Effects (IEs) associated with a phenomenon. Many studies determine ε using a Rayleigh-plot, which needs multiple data points. More recent articles describe an alternative method using the Rayleigh equation that allows the determination of ε using only one experimental point, but this method is often subject to controversy. However, a calculation method using two points (one experimental point and one at t0) should lead to the same results because the calculation is derived from the Rayleigh equation. But, it is frequently asked "what is the valid domain of use of this two point calculation?" The primary aim of the present work is a systematic comparison of results obtained with these two methodologies and the determination of the conditions required for the valid calculation of ε. In order to evaluate the efficiency of the two approaches, the expanded uncertainty (U) associated with determining ε has been calculated using experimental data from three published articles. The second objective of the present work is to describe how to determine the expanded uncertainty (U) associated with determining ε. Comparative methodologies using both Rayleigh-plot and two point calculation are detailed and it is clearly demonstrated that calculation of ε using a single data point can give the same result as a Rayleigh-plot provided one strict condition is respected: that the experimental value is measured at a small fraction of unreacted substrate (f < 30%). This study will help stable isotope users to present their results in a more rigorous expression: ε ± U and therefore to define better the significance of an experimental results prior interpretation. Capsule: Enrichment factor can be determined through two different methods and the calculation of associated expanded uncertainty allows checking its significance.

Position-Specific 13C Fractionation during Liquid-Vapor Transition Correlated to the Strength of Intermolecular Interaction in the Liquid Phase.

The relationship between the strength of the intermolecular interaction in liquid and the position-specific 13C fractionation observed during distillation was investigated. A range of molecules showing different intermolecular interactions in terms of mode and intensity were incorporated in the study. Although it had previously been suggested that during evaporation the diffusive 13C isotope effect in the thin liquid layer interfaced with vapor is not position-specific, herein we show that this is not the case. In particular, the position-specific effect was demonstrated for a series of alcohols. Our hypothesis is that intermolecular interactions in the liquid phase are the source of position-specific 13C fractionation observed on the molecule. A clear trend is observed between the 13C isotope effect of the carbon bearing the heteroatom of chemical function and the relative permittivity, the solvent hydrogen bond acidity, and the solvent hydrogen bond basicity, while only a weak trend was observed when using the 13C content of the whole molecule. Furthermore, two families of products appeared when using the hydrogen bond acidity parameter for the correlation by distinguishing H-acceptor and H-donor molecules from those H-acceptors only. This strongly reinforces the hypothesis of an important role of the 13C positioned close to the interaction center.

Simulating Stable Isotope Ratios in Plumes of Groundwater Pollutants with BIOSCREEN-AT-ISO.

BIOSCREEN is a well-known simple tool for evaluating the transport of dissolved contaminants in groundwater, ideal for rapid screening and teaching. This work extends the BIOSCREEN model for the calculation of stable isotope ratios in contaminants. A three-dimensional exact solution of the reactive transport from a patch source, accounting for fractionation by first-order decay and/or sorption, is used. The results match those from a previously published isotope model but are much simpler to obtain. Two different isotopes may be computed, and dual isotope plots can be viewed. The dual isotope assessment is a rapidly emerging new approach for identifying process mechanisms in aquifers. Furthermore, deviations of isotope ratios at specific reactive positions with respect to "bulk" ratios in the whole compound can be simulated. This model is named BIOSCREEN-AT-ISO and will be downloadable from the journal homepage.

Simulating stable carbon and chlorine isotope ratios in dissolved chlorinated groundwater pollutants with BIOCHLOR-ISO.

BIOCHLOR is a well-known simple tool for evaluating the transport of dissolved chlorinated solvents in groundwater, ideal for rapid screening and teaching. This work extends the BIOCHLOR model for the calculation of stable isotope ratios of carbon and chlorine isotopes in chloroethenes. An exact solution for the three-dimensional reactive transport of a chain of degrading compounds including sorption is provided in a spreadsheet and applied for modeling the transport of individual isotopes 12C, 13C, 35Cl, 37Cl from a constant source. The model can consider secondary isotope effects that can occur in the breaking of CCl bonds. The model is correctly reproducing results for δ13C and δ37Cl modeled by a previously published 1-D numerical model without secondary isotope effects, and is also reproducing results from a microcosm experiment with secondary chlorine isotope effects. Two applications of the model using field data from literature are further given and discussed. The new BIOCHLOR-ISO model is distributed as a spreadsheet (MS EXCEL) along with this publication.